arXiv:1507.03771vl [quant-ph] 14Jul2015 


Fast bias inversion of a double well without residual particle excitation 
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We design fast bias inversions of an asymmetric double well so that the lowest states in each well 
remain so and free from residual motional excitation. This cannot be done adiabatically, and a 
sudden bias switch produces in general motional excitation. The residual excitation is suppressed 
by complementing a predetermined fast bias change with a linear ramp whose time-dependent slope 
compensates for the displacement of the wells. The process, combined with vibrational multiplexing 
and demultiplexing, can produce vibrational state inversion without exciting internal states, just by 
deforming the trap. 

PACS numbers: 32.80.Qk, 03.Be, 37.10.Gh, 37.10.Vz 


I. INTRODUCTION 

In a recent paper [l| a protocol to realize fast vi¬ 
brational state multiplexing/demultiplexing of ultra cold 
atoms was introduced. By a properly designed time- 
dependent potential deformation between a harmonic 
trap and a biased double well, the states of a single 
atom in a harmonic trap can be dynamically mapped into 
states localized at each well (vibrational demultiplexing, 
see the left arrow in Fig. 1), or vice versa (multiplex¬ 
ing, see the right arrow in Fig. 1), faster than adiabat¬ 
ically and without residual excitation at the final time. 
It was suggested that these processes may be combined 
with a bias inversion to produce state inversions, from 
the ground to the hrst excited state of the harmonic trap 
and vice versa, based on trap deformations only, see Fig. 
[U or to produce vibrationally excited Fock states from an 
initial ground state [I| . These are basic operations to im¬ 
plement quantum information processing. Thus the pos¬ 
sibility to perform them without exciting internal atomic 
states as an intermediate step is of much interest, as deco¬ 
herence induced by decay would be suppressed; moreover 
the systems amenable of manipulation would not need 
to have an isolated two-level structure. The objective 
of this paper is to design fast controlled bias inversions 
so that the lowest states in each well remain so without 
residual excitation. Unlike multiplexing, however, there 
is no truly adiabatic slow process that achieves this state 
transformation. In the bias inversion depicted within the 
central frame of Fig. [U for example, a slow bias change 
would preserve the state ordering so that the states rep¬ 
resented in the third potential configuration would be 
interchanged (i.e. the grey state in the right well, and 
the white one in the left well). Nevertheless, in the limit 
in which the two wells are effectively independent, which 
in practice means, for times shorter than the tunneling 
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FIG. 1: Schematic representation of demultiplexing (left ar¬ 
row), bias inversion (framed in dashed line, central arrow) 
and multiplexing (right arrow). 


time among the wells, the intended state transition might 
indeed be done slowly enough to be considered adiabatic. 
If we approximate each “isolated” well by a harmonic os¬ 
cillator, the intended transformation amounts to a “hor¬ 
izontal” displacement along the inter well-axis together 
with a rising/lowering of the energy of the wells. The 
latter effects, however, do not affect the intrawell dy¬ 
namics, so we may focus on the displacement. In other 
words, within the stated approximations the bias inver¬ 
sion amounts to the transport of a particle in a harmonic 
potential. Thus, to achieve a fast transition without 
residual excitation we may use shortcuts to adiabatic- 
ity designed to perform fast transport 0 . Specifically we 
shall use a compensating-force approach [3,I3, equivalent 
to the fast-forward scaling technique @ , based on adding 
to the potential a linear ramp with time-dependent slope 
to compensate for the effect of the trap motion in the 
non-inertial frame of the trap. We shall compare this ap¬ 
proach with a sudden bias switch, a fast quasi-adiabatic 
(FAQUAD) approach or a smooth polynomial con¬ 
nection without compensation. In Sec. El we introduce 
the compensating-force approach and Sec. IIIII describes 
the alternative methods. In Sec. IIVI numerical examples 
are presented with parameters appropriate for trapped 
ions in multisegmented, linear Paul traps, and for neu¬ 
tral atoms in optical traps. Finally, in Sec. |V]we discuss 
the results and open questions. 
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II. COMPENSATING-FORCE APPROACH 


III. ALTERNATIVE METHODS 


If the double well potential with nearly independent 
wells is subjected to a bias inversion such that the trap 
frequencies of each well are essentially equal and constant 
throughout, and the trajectories of the well minima move 
in parallel, the process may be described by a parallel 
transport of two rigid harmonic oscillators, one for each 
well. The Hamiltonian potential near the minima may 
be approximated as 


Vq{x- xo) = ]^mnl{x-xof, (1) 

where Hq is the angular frequency and xq = xo(t) is the 
common notation for either of the two minima.^ When 
needed we may distinguish the minima as a;o,±, with 
> xo,-. The Hamiltonian Hq = /2m + Vo(x — xq) 

has eigenenergies En = {n + ^) hQo, and well known nor¬ 
malized eigenstates 4>n{x — xq), proportional to Hermite 
polynomials Q- 

Adding to the Hamiltonian a linear term with an ap¬ 
propriate time-dependent slope the non inertial effect of 
the motion of the well will be compensated in the trap 
frame ii. To define the trap frame consider the follow¬ 
ing position and momentum displacement unitary oper¬ 
ator 


^ _ ^ipxo(J.)/h^-imxo{t)xlh 

where the overdot represents a time derivative. Starting 
from the Schrodinger equation 

ihdt\tp) = Hq\iP), (3) 

the transformed wave function |$) = U\tjj) obeys 

ihdt\<^) = UHoU^‘i>) + in{dtU)U'<\^) = ( 4 ) 

where the interaction picture (trap frame) Hamiltonian 
is 


H/ = 


2m 


+ Vo{x) + mxxo + TTwii 


0i 


( 5 ) 


and Vq{x) = ^m^l^x^. The term mx/^/2 only depends on 
time, it does not affect the dynamics and can be ignored. 
To compensate the motion of the trap, we add —mxxo 
to Hf). This produces —m(x + xo)io in the trap frame 
and the resulting potential in that frame is reduced to 
Vo(x), again neglecting purely time-dependent functions. 
Vo{x) does not depend on time, so any stationary state 
in this trap frame will remain stationary, an excitations 
are avoided. 


^ We disregard purely time-dependent functions in each well. Dif¬ 
ferential phases among the wells depending on these functions 
can be ignored since the traps are assumed to be independent. 


In this section we present three simple alternative ap¬ 
proaches to perform the bias inversion. 

Sudden approach. In the sudden approach the poten¬ 
tial is changed abruptly from the initial to the final con¬ 
figuration, but the state of the system remains unchanged 
immediately after the potential change (in general it will 
evolve afterwards). If the target state is i/tar the result¬ 
ing fidelity is 


Fs = \{tp{0)\tj}tar)\- (6) 

Fast quasi-adiabatic approach (FAQUAD). A quasi- 
adiabatic method to speed up adiabatic processes when 
there is one control parameter X{t) is based on distribut¬ 
ing the adiabaticity parameter homogeneously in time 
0. For instantaneous levels 0 and 1 this means 


h 


Eq — El 


( 7 ) 


where the instantaneous eigenstates (po, (/i eigenen¬ 
ergies Eg, El depend on time through their dependence 
on A, and c is constant. By the chain rule this becomes 
a first order differential equation for A(t), and c is set so 
that the boundary conditions for A(t) at initial, t = 0, 
and final time tf are satisfied. In the transport of a par¬ 
ticle with a harmonic oscillator of angular frequency Hq 
centered at xo{t) we set X{t) = xo{t). Using the energies 
and eigenstates of the first two levels of the harmonic 
oscillator in Eq. ©, the FAQUAD condition becomes 
simply 


mxo(t) 

v^2ftmDo 


( 8 ) 


The solution is the linear connection X(i{t) = a;o(0) -I- 
[xo(tf) — xo(0)]t/tf. The minimal tj for which a zero of 
excitation energy appears is @,[3 27r/f2o. 

Polynomial connection without compensation. The fi¬ 
nal and initial values of the control parameter may as 
well be smoothly connected without applying any com¬ 
pensation, for example using a fifth order polynomial that 
assures the vanishing of first and second derivatives of the 
parameter at the boundary times. 


IV. EXAMPLES 

In the following examples, the potentials and parame¬ 
ters are adapted for a trapped ion in a multisegmented 
Paul trap, and for a neutral atom in a dipole trap. 


A. Trapped ions 

For a trapped ion we consider a simple double well 
potential of the form 

V{x, t) = fdx^ ax^ 'yx, 


( 9 ) 
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with a(t) < 0 and j3{t) > 0 iSIl- a and P are assumed 
to be constant and 7 = 7 (t) time dependent. The bias 
inversion implies a change of sign of 7 (t) from 70 > 0 to 
—7o- 

From ^ = 0 the condition for the extrema becomes 


+ 2ax + 7 = 0 . 

It is useful to define 


( 10 ) 


( 11 ) 


-3B 2A^ - 9AB + 27C 

Q — n ’ ^ — ~A • 


9 


54 


When B? < there are two minima and one maximum. 
With a < 0 and /3 > 0, this is satisfied for 


I 7 I < 


9 


3/2 


The trajectories of the minima are 

' + (l±l)7r\ A 


Xo 


,± = -2V^ 


COS 


3 ’ 


(13) 


(14) 


R 

taken to be positive. Up to second order in 7 they are 



1 

1 a 7 37 ^-\/—a/3 

(15) 

Xo- - 

“ ■72\ 

/ P Aa lQy/2a^ ’ 


I r 

a 1 7 1 ij'^yf-aP 

(16) 


“ 71V’ 

/3 4a 1642a3 


^, 0 < 0 < TT and the roots are 


where 9 = arccos 


The distance travelled by each well is, when (ITTl) is ful¬ 
filled, d = 70 /(2a), see Eqs. (fT^ and (ITCl) . and the 
effective frequency at each minimum 


Wo = 


1 UW 


-\dx^Jx=xo' 
For Eq. © the effective frequencies are 


( 20 ) 


Wo,: 


,± = \ —\ a+-P <A+6y^cos 


lfe+{l± l)7r' 

3 




m 2\p2 V a^m 


-7- 


( 21 ) 


Hence, comparing the two terms, the condition for the 
frequencies to be essentially constant wo,=p ~ Uq = 
2 \/~ m again the inequality in Eq. ( 113 - 
In the regime where the inequality (11711 holds, we can 
apply the compensating force approach to implement a 
fast bias inversion. Since the compensating term —mxio 
is equal for both harmonic traps, we add it to V in Eq. 
(|9]), and the resulting Hamiltonian JI is 


H = 


2m 


+ Px'^ + ax^ + (7 ~ mxo)x. 


( 22 ) 


Note that the compensation amounts to changing the 
time dependence of the slope of the linear term from 
the reference process defined by 7 (t) to 7 e//(i) = 'lit) — 
mxoit) = 7 (t) — m 7 (t)/( 4 a). To set 7 (t) we design a 
connection between the initial and final configurations. 
Eirst note the boundary conditions 


7 ( 0 ) = 70 > 0, litf) = - 70 , (23) 


The quadratic term 7 is negligible with respect to the 
linear term when 


|7| < 



(17) 


which implies that the trajectories for the minima move 
in parallel. Note that this inequality automatically im¬ 
plies the one in Eq. m- Neglecting the quadratic term, 
the two minima are given by a;o,± = 

The distance between the minima is given by 


which we complement by 

7 ( 4 ) = 0, 7 ( 4 ) = 0, (24) 

4 ~ ^'>tfj 

so that xoitb) = io(4) = 0. This implies that W(tb) = 
^ipxoiH)/h Uitb) = 0. Therefore the Hamiltonians 
and the wave functions in interaction and Schrbdinger 
pictures transform into each other by a simple coordinate 
displacement. At intermediate times, we interpolate the 
function as jit) = where the coefficients are 

found by solving Eqs. ((^51) and (|24l) . Finally, 


D 



sm 


:i7T + 29) 



W-aP 2 
8V2a^ ^ 


(18) 


We may also compute the energy bias between the two 
wells as 


(19) 


lit) = 7 ( 0 ) + 10 ( 7 ( 4 )-7(0)]s^ 

- 15(7(4)-7(0)]s^ +6(7(4)-7(0 )]s4 ( 25 ) 

where s = t/tf. This function and examples of Jeff are 
shown in Fig. [2] 

In order to compare the robustness of the compensat¬ 
ing force method against the alternative ones we consider 
a ®Be+ ion in the double well with the realistic param¬ 
eters a = —4.7 pN/m and P = 5.2 mN/m^, similar to 


5 = jD. 
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FIG. 2: (Color online) 7 versus s = t/tf for the polynomial 
in Eq. 11251) (solid black line) and FAQUAD (short-dashed 
red line). 70 = 86.4 zN, 7 (t/) = — 70 , a = —4.7 pN/m, 
and /I = 5.2 mN/m®. Also shown are the different effective 
slopes adding a compensation to the polynomial, 7e//(t) = 
7 (t) — m 7 (t)/( 4 a), for the mass of ®Be’'" and times tf = 0.07 
fis (long-dashed blue line); tf = 0.1 fj,3 (green dots); and 
tf = 0.3 fis (magenta squares). 

those in [l^. For a moderate initial bias compared to 
the vibrational quanta, such as 


the fidelity provided by the sudden approach is one for all 
practical purposes so we can change the bias abruptly and 
reach the target state. The displacement of the trap d 
may be compared with the oscillator characteristic length 
tto = ^Jh/mQ,Q. Their ratio is 


For the Paul trap R 0.00065, which explains the high 
fidelity of the sudden approach for a moderate bias in¬ 
version of the ion. At these bias values there is really 
no need to apply a more sophisticated protocol than the 
sudden switch. 

Henceforth, we assume a much larger 70 , but still sat¬ 
isfying the condition (1171) . In particular, for an initial 
bias of 1000 flo?i (corresponding to 70 = 86.4 zN), the 
ratio becomes R ~ 0.65. The maximum variation of the 
difference between the trajectories of the minima is 3 
pm, about three orders of magnitude less than the dis¬ 
placement of each minimum (9.2 nm), so the minima 
follow parallel trajectories. Furthermore, the maximum 
variation of the frequencies in Eq. (I^Tl) with respect to 
Ho = 27r X 5.6 MHz is 27r x 3.7 kHz, so the effective 
frequency is nearly constant. 

Figure [3] demonstrates the effect of the compensating- 
force approach. Starting from the ground state of the 
lower (left) well, the final evolved state following the 
shortcut with compensation stays as the “ground state” 
of the left well. This is actually defined as the lowest 
state of the double well system predominantly located 



- 21.2 - 21.1 21.1 21.2 
X (pm) X (pm) 


FIG. 3: (Color online) Left: Ground state of the left well 
at t = 0 (long-dashed magenta line) and ai t = tf (magenta 
triangles), and final state with the compensating force applied 
on the double well (solid blue line). Right: Ground state of 
the right well: at t = 0 (short-dashed red line) and at t = tf 
(red dots) and final state with the compensating force applied 
(solid black line), tf = 4 ns and other parameters as in Fig. 
El for ®Be+. 



FIG. 4: (Color online), (a) Fidelity \{4>Litf)\'ilj{tf))\, where 
\4‘L{tf)) is the lowest state located in the left well in the final 
time configuration, and |'i/'(t/)) is the evolved state following 
the shortcut at final time, (b) Final excitation energy for 
the process on the left well using compensating-force (blue 
dots), fifth degree polynomial in Eq. (1251) (solid black line), 
and FAQUAD (short-dashed red line). The parameters are 
for ®Be+ as in Fig. E) 


on the left. There is a similar process for the right well. 
The final states represented are obtained by solving the 
Schrodinger equation with the full Hamiltonian (|2^ . Fig. 
0] is for the process in the left well. It compares the fi¬ 
delity at final time and the excitation energy, defined as 
Eex = E{tf)—EQ{tf) where E{tf) is the final energy after 
the quantum evolution following the shortcut and E^itf) 
is the ground state final energy of the upper harmonic 
well, using the polynomial (1251) for 7 with and without 
compensation, as well as the results of the FAQUAD ap¬ 
proach. The fidelity without compensation tends to the 
fidelity of the sudden approach (0.89) for very short final 
times. The method with compensation clearly outper¬ 
forms the others. In principle a fundamental limitation 
of the approach is due to the fact that the inequality dni, 
that guarantees parallel motion and stable frequencies of 
the wells, should as well be satisfied by 7 e//, but, at very 
short times, the dominant term of 7 e// ~ —mq/da may 
be too large. To estimate this short time limit we com¬ 
bine the mean-value theorem inequality for the maximum 
0, \l\max > with Eq. dni) for 7 e// to find the 
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condition 


The factor on the right hand side is 10s for this ex¬ 
ample, see Fig. |4l which is several orders of magnitude 
smaller than 27r/r2o and does not affect the fidelity in the 
scale of times shown. 


B. Neutral atoms 


The potential 


Vnaix, t) = -I- Vq COS^ 


7r(x — Ax) 
di 


(29) 


forms also a double well. It was implemented in [T^ with 
optical dipole potentials, combining a harmonic confine¬ 
ment due to a crossed beam dipole trap with a periodic 
light shift potential provided by the interference pattern 
of two mutually coherent laser beams, uj is the angular 
frequency of the dipole trap about the waist position, Vq 
the amplitude. Ax the displacement of the optical lat¬ 
tice relative to the center of the harmonic well, and di 
is the lattice constant. (Double wells with a controllable 
bias may be also realized by two optical lattices of differ¬ 
ent periodicity with controllable intensities and relative 
phase 0 .) Here, the bias inversion implies a change of 
sign of Ax{t) from Accq > 0 to —Axq. 

To check if the conditions to apply the compensating 
force approach hold here as well, an analysis similar to 
the one in the previous example is now performed. We 
approximate the potential around each minimum, 
up to fourth order. From = 0 we get a cubic 

equation for each minimum. The positions of the minima 
are thus given by 


^0,± 



- 27r 

3 


H(±) 

3 


(30) 


where 


Q 

^(±) 

0(±) 


2d^TT^VQ + dfmuj'^ 
dTT^Vb ’ 

-^{2Ax ±di), 

—3di{2Ax ± dAm'K'^\/V quj^ 

pf)0 _1_ i _ 

2(27r2Fo + d2mw2)3/2 


Up to a second order in Acc, 


a;o,± ~ ±a-I-foAa; ± cAa;^, (32) 


where the coefficients are known explicitly but too 
lengthy to be displayed here. Whenever the quadratic 
term is negligible with respect to the linear term 


(cAa;^ << bAx), we can approximate Xo,± = ±a + foAa; 
(parallel motion). The distance between the minima is 
given by 


D 


1 

3 


+ 6v^ 


— COS 




2a -|- 2cAx^. 


(33) 


Moreover, a;o,± ~ f^gAx, again with known but lengthy 
coefficients g and /. As long as gAx << / we may set 
wo,± ~ Ho = /. 

For realistic parameters the conditions for parallel 
transport of the minima and constant frequency are in¬ 
deed satisfied. We consider a ®^Rb atom in the trap and 
set the parameters in [ij after the demultiplexing process, 
namely, di = 5.18 fim, to = 2 tt x 59.4 Hz, and Vb//i = 1.4 
kHz; the time-dependent displacement Ax = Ax(t) is 
the control parameter to perform the bias inversion, such 
that 


Aa;(0) = Axo, ^x(tf) = —Axq, (34) 

with Axq = 200 nm. We also impose 

Aa;(0) = 0, Aa;(0) = 0, 

Ax(tf) = 0, Ax(tf) = 0, (35) 


to achieve similar conditions in the derivatives of the 
minima xq- At intermediate times, we interpolate the 
function as Ax(t) = where the coefficients 

are found by solving Eqs. o and (ESD. Consequently, 
the connection between the initial and final Hamiltonians 
is given by the same polynomial in Eq. (j25p changing 
7 (f) —>■ Ax{t). The double wells are much deeper and 
tight for trapped ions than for neutral atoms, compare 
an intrawell angular frequency Hq of 27r x 5.6 MHz for the 
ions versus 27r x 0.35 kHz for the optical trap. Therefore, 
in this case, for a moderate initial bias compared to the 
vibrational quanta, the ratio between the displacement 
of the trap d and the oscillator characteristic length oq is 
i? ft! 0.67. 

With the parameters given at time f = 0, the sepa¬ 
ration of the minima is D = 5 ^m, the bias between 
minima S = 2.02 x 10“32 and an the effective angular 
frequency Hq = 27r x 0.35 kHz, whereas the maximum 
variation of the frequencies along the process is approx¬ 
imately 27r X 0.2 Hz. Furthermore, the maximum devia¬ 
tion from D of the minima separation is 1.8 nm, whereas 
the displacement of each minimum is about 0.4 /rm. In 
summary the minima do move in parallel with constant, 
equal frequencies for practical purposes. 

To accelerate the bias inversion we add the compen¬ 
sating term to V in Eq. (EH), 


H = -^-t-mw^x^-l-Vo cos^ 
2m 2 


t:(x — Ax) 
di 


—mxxQ. (36) 


Figure [5] shows the evolution of the densities. Focusing 
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l¥(x.t)|^ 6 0.8 



-4 -2 0 2 4 

X (|im) 


FIG. 5: (Color online) Evolution of the wave function densi¬ 
ties following the shortcut in Eq. (IMll for states in left and 
right wells. The parameters are for ®^Rb: di = 5.18 pm, 
uj = 59.4 X 27r Hz, Vojh = 1.4 kHz, Axo = 200 nm, and 
tf — 63 ps. 



FIG. 6: (Golor online), (a) Fidelity |(<pi(t/)|V'(i/))li where 
\<pi{tf)) is the lowest state predominantly of the left well at fi¬ 
nal time (the first excited state of the double well) and |V'(l/)) 
is the evolved state following the shortcut at final time, (b) 
Final excitation energy. Gompensating-force approach (blue 
dots), fith degree polynomial in Eq. (I25II with the change 
7 (t) —>■ Ax{t) without compensation (solid black line), and 
FAQUAD approach (short-dashed red line). The parame¬ 
ters are chosen for ®^Rb: di = 5.18 pm, cj = 59.4 x 27r Hz, 
Vo/h = 1.4 kHz, and Aiq = 200 nm. 


on the left well. Fig. [5] (a) demonstrates that the fidelity 
is exactly one (blue dots) with the compensating force. 
However, using for the inversion the fifth degree polyno¬ 
mial in Eq. (ES|) (with the change 7 ( 1 ) —?► Ax{t)) without 
compensation the fidelity at short final times decreases 
until the value of the sudden approach, 0.9. Furthermore, 
the excitation (residual) energy is approximately 
zero for the shortcut protocol, compared to the excitation 
without compensation in Fig. in](b). In practice, the min¬ 
imal process times tf are not limited by the method per se 
but by the technical capabilities to implement the maxi¬ 
mal compensating force. This force depends on the maxi¬ 
mal acceleration of the well, whose lower bound is known 
to be Umax = Q. To implement the compensation 

with a magnetic field gradient G, the gradient should be 
of the order of G ~ inamax /in an amount of time tf 
(ps is Bohr’s magneton). For rubidium atoms polarized 


in the magnetic level F = mp = 2 and the transport pa¬ 
rameters considered above, this requires a magnetic field 
gradient on the order of 40 T/m shaped on a time in¬ 
terval tf = 63 ps. This is definitely challenging from an 
experimental point of view. Alternatively, one can use 
the dipole force of an out of axis Gaussian laser beam. If 
the double well is placed on the edge at a distance of w /2 
where w is the waist and if ap denotes the polarizabil¬ 
ity, the local potential slope experienced by the atoms 
is on the order of apP/w^ where P is the power of the 
beam. The compensation requires that P/w^ = mjap. 
For instance, with an out-of-resonance beam at a wave 
length of 1 pm, the polarizability of rubidium-87 atoms 
is ap ~ 1.3 X 10“^® m^s, and the compensation can be 
performed using a IW laser with a waist of 20 pm. 


V. DISCUSSION 


In this work we have proposed a method to invert 
the bias of a double-well potential, in the regime of 
independent wells, to keep the final states motionally 
unexcited within the same original well. The method 
treats the bias inversion as a rigid transport of the wells, 
which is justified for realistic parameters, and applies a 
“compensating-force” to cancel the excitations. Exam¬ 
ples have been worked out for ions or neutral atoms and 
comparisons have been provided with a sudden approach, 
a fast quasi-adiabatic (FAQUAD) approach, or a smooth 
polynomial connection of initial and final bias without 
compensation. The compensating-force method clearly 
outperforms the others and gives ideally perfect fidelities, 
at least under ideal conditions, up to very small times 
compared e.g. with the time 27r/Do where FAQUAD pro¬ 
vides a first zero of excitation. 

One more advantage is the flexibility, as the reference 
process used to design the corresponding compensation 
(we have used a polynomial for simplicity) may be cho¬ 
sen among a broad family of functions satisfying Eqs. 
(1^ [23|). As in other shortcut approaches, this flexibil¬ 
ity may be used to enhance robustness versus noise and 
perturbations [IMl- 

The bias inversion put forward here and the multiplex¬ 
ing and demultiplexing protocols developed in [I| provide 
the necessary toolkit to perform vibrational state inver¬ 
sions [H, [ 1 ^, or Fock state preparations [l|. Applica¬ 
tions in optical waveguide design are also feasible [^ . 
As well, the fast bias inversion is directly applicable to 
Bose-Einstein condensates [2l|, (2^ . 

Generalizations for conditions in which rigid transport 
does not hold are also possible using invariant theory Q , 
which allows for finding processes without final excitation 
where both the frequency and position of the well depend 
on time (^ . 
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